An initialization subroutine for DAEs solvers: DAEIS
نویسنده
چکیده
Numerical difficulty arises when DAEs are solved with inconsistent initial values of dependent variables. It can cause the solution failures of many popular DAEs solvers. An initialization subroutine, DAEIS (Differential Algebraic Equations Initialization Subroutine), has been developed to handle this issue. In DAEIS, consistent initial values of dependent variables for index-1 DAEs are obtained through solving a set of nonlinear equations based on the initialization settings. To achieve good efficiency and robustness of the initialization computation, a nonlinear equation solver, GNES, has been specially built for DAEIS. The usage of DAEIS is demonstrated through handling inconsistent initializations for combined continuous/discrete simulations. © 2001 Elsevier Science Ltd. All rights reserved.
منابع مشابه
Extending explicit and linearly implicit ODE solvers for index-1 DAEs
Nonlinear differential-algebraic equations (DAE) are typically solved using implicit stiff solvers based on backward difference formula or RADAU formula, requiring a Newton-Raphson approach for the nonlinear equations or using Rosenbrock methods specifically designed for DAEs. Consistent initial conditions are essential for determining numeric solutions for systems of DAEs. Very few systems of ...
متن کاملIndex Reduction and Discontinuity Handling using Substitute Equations
Several techniques exist for index reduction and consistent initialization of higher index DAEs. Many such techniques change the original set of equations by differentiation, substitution, and/or introduction of new variables. This paper introduces substitute equations as a new language element. By means of a substitute equation, the value of a continuous variable or its time derivative is spec...
متن کاملAdjoint Sensitivity Analysis for Differential-Algebraic Equations: The Adjoint DAE System and Its Numerical Solution
An adjoint sensitivity method is presented for parameter-dependent differentialalgebraic equation systems (DAEs). The adjoint system is derived, along with conditions for its consistent initialization, for DAEs of index up to two (Hessenberg). For stable linear DAEs, stability of the adjoint system (for semi-explicit DAEs) or of an augmented adjoint system (for fully implicit DAEs) is shown. In...
متن کاملConsistent Initial Conditions for Unstructured Higher Index DAEs: A Computational Study
Differential algebraic equations (DAEs) are implicit systems of ordinary differential equations, F (x′, x, t) = 0, for which the Jacobian Fx′ is always singular. DAEs arise in many applications. Significant progress has been made in developing numerical methods for solving DAEs. Determination of consistent initial conditions remains a difficult problem especially for large higher index DAEs. Th...
متن کاملTowards Integrating Hybrid DAEs with a High-Index DAE Solver
J.D. Pryce and N.S. Nedialkov have developed a Taylor series method and a C++ package, Daets, for solving numerically an initial-value problem differential-algebraic equation (DAE) that can be of high index, high order, nonlinear, and fully implicit. Numerical results have shown this method to be efficient and very accurate, and particularly suitable for problems that are of too high an index f...
متن کامل